Searching for the Sources Responsible for 
Cosmic Reionization: Probing the Redshift 
Range 7<z<10 and Beyond 



Daniel P. Stark 1 k Richard S. Ellis 1 

l Dept. of Astronomy, California Institute of Technology, Pasadena CA 91125 

E-mail: dps@astro.caltech.edu 



Abstract 

We review recent observations that suggest that the global star formation rate 
density of UV-bright galaxies is declining monotonically with redshift over 3< z <7 
and illustrate the challenges that a continuation of this decline poses in explaining 
the assembled stellar mass in several z~6 galaxies deduced from recent Spitzer data. 
A plausible conclusion is a vigorous period of yet earlier star formation. Prior to 
JWST and TMT, strong lensing offers a unique probe of the extent of this earlier 
activity. We discuss the first results of a blind spectroscopic survey of 10 lensing 
clusters for 8.5< z <10 Lya emitters using NIRSPEC on Keck. We demonstrate 
this survey is achieving an (unlensed) sensitivity equivalent to a star formation rate 
of 0.1 M & yr _1 . A companion survey, now beginning with HST and Spitzer, will 
target lensed z and J-band dropouts and probe a ~1 arcmin 2 region 1 magnitude 
deeper than the UDF/NICMOS observations. In combination, both surveys will 
provide the first constraints on the contribution of early, low luminosity, sources to 
cosmic reionization. 
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1 Introduction 



The reionization of neutral hydrogen which rendered the Universe transparent 
to UV photons and thereby terminated the so-called 'dark ages' was a land- 
mark event in cosmic history. Determining the nature, abundance and redshift 
distribution of the sources responsible for this important transition represents 
the current frontier in studies of the first galactic systems. 

Recent developments strongly motivate a systematic search for star-forming 
sources beyond z~7. The WMAP temperature-polarization cross correlation 
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signal on large angular scales (Jlj) implies scattering of microwave photons by 
free electrons from ionizing sources at z ~10-12 (J2J). The abundance of lumi- 
nous UV-emitting i— band (ji)), z— band, and J— band (jl),(0) 'drop-outs' in the 
Hubble Ultra Deep Field at z~6, 7 and 10 respectively, suggests a surprisingly 
continuous decline in the comoving density of emerging UV photons to z=10. 
The luminosity density of these UV-bright sources seems insufficient to reion- 
ize the universe unless non-standard initial mass functions (0) or very steep 
faint end-slopes for the luminosity function (J7j) are invoked. More likely is the 
possibility of an intense earlier period of star formation. Meanwhile, Spitzer 
observations have uncovered a population of very massive (M^ ar ~ 10 10 - 
10 11 M ) galaxies with well-established stellar populations flih. (|2oT). 

(2J). Given the decline in the star formation rate density of UV-emitters, the 
presence of this massive population requires explanation. While extinction and 
cosmic variance likely play a role in the apparent discrepancy, an earlier period 
of vigorous star formation may have contributed significantly to buildup these 
massive galaxies as well as producing a substantial fraction of the photons 
required for reionization. 



Finding this early population is one of the primary goals of observational cos- 
mology. With current facilities, traditional search methods are poorly equipped 
to search beyond z~7. Even with HST and 8-10 meter class telescopes, con- 
tinuum ('drop-out') and narrow-band imaging Lyman a searches only probe 
the most luminous, rare, sources at these redshifts; the contribution to the 
reionization flux from these sources may well be small. The study of lower lu- 
minosity sources rendered visible by the strong gravitational lensing induced 
by clusters with well-constrained mass models can greatly assist in this work. 
A typical lensing cluster magnifies sources by x5-10 over 1 arcmin 2 for sources 
at z >7. Moreover, on the critical line, faint signals are boosted >x30 and 
spectroscopy and more detailed studies are often feasible at otherwise impos- 
sible (unlensed) limits. For several objects recently located via this technique 
W)i ficft- stellar continuum slopes, Lya profiles and star formation rates 
have been determined for sources whose intrinsic flux is close to the UDF 
limit, Iab —30. Perhaps most remarkably, such sources have been magnified 
so as to become accessible by Spitzer fjllh yielding luminosity-weighted ages. 



Until the era of TMT/JWST, scanning the critical lines of massive clusters 
provides possibly the only reliable means of estimating the abundance of low 
luminosity star forming sources in the 7< z <10 interval with minimal con- 
tamination. In these proceedings, we discuss two such surveys: in §4, we report 
on the progress of a Keck spectroscopic survey along cluster critical lines for 
Lya emitters at z>8, and in §5 we present a new 110 orbit Hubble imaging 
survey for lensed low-luminosity galaxies at z~7-10. 
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2 Evidence for a Declining Star Formation Rate 



Optical and near-infrared imaging with the Hubble Space Telescope have ten- 
tatively provided the first limits on the visible star formation activity in the 
era 6 < z < 10. Galaxies are identified at these redshifts by extending the 
Lyman-break technique fl3). (flih to z>6, and their star formation rates are 
estimated by measuring their rest-frame ultraviolet emission, which is domi- 
nated by young massive stars. 

The abundance of z~6 galaxies (identified as i'-dropouts) has been calculated 
in the 400 arcmin 2 GOODS-North and GOODS-South fields down to z AB <25.6 
(15), in ACS parallel observations covering 10 arcmin 2 down to z^ B <28.3 
(17); i n ACS GTO fields spanning 46 arcmin 2 reaching down to zab<27.3 
( la), in the ACS parallels to the NICMOS UDF observations which span 21 
arcmin 2 down to z^s<28.1, and in the 11.3 arcmin 2 UDF down to zab < 28.5 
(0) and down to zab < 29.2 (0). In addition, the exceptional depth of the HST 
NICMOS images of the UDF allows galaxies at z~7 and z~10 to be identified 
as z-drops and J-drops, respectively (jl) , (jfj). 

The general consensus emerging from these studies is that the apparent star 
formation rate density, uncorrected for extinction (integrated down to L[/y>0.1 
L* =3 ) is ~3-6 times less at z~6 than at z~3 and continues to decline out to 
z~10 (integrated down to L[/y>0.3L* =3 ) (Fig 1). The declining star forma- 
tion rate poses a challenge for models suggesting star-forming galaxies at z~6 
reionized the Universe. The low-abundance of high-z dropouts can be recon- 
ciled with such models if the galaxies have low metallicity or a non-standard 
initial mass function (0), or if the UV luminosity function has a very steep 
faint-end slope (J7J), such that the bulk of ionizing photons are released by 
intrinsically faint galaxies that lie below the 5a sensitivity limit of the UDF. 
Alternatively, the bulk of reionization may occur at yet earlier times, requir- 
ing the presence of a significant abundance of objects at z^>6, or conceivably 
reionization was not accomplished by young star-forming galaxies. 



3 The Crucial Role of Spitzer 



Detecting light from the most distant galaxies known in the universe has been 
one of the most remarkable and important feats of the IRAC instrument on- 
board the 85 cm Spitzer Space Telescope. The IRAC filters at 3.6-8.0 /j,m probe 
the rest-frame optical and near-infrared at z>6, providing the first glimpse of 
light from established stellar populations of galaxies at such high redshift. 
Combining these data with deep broadband optical photometry from HST 
and 8-10 meter class ground based telescopes, spectral energy distributions 
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Fig. 1. The co-moving volume density of luminous UV emitting galaxies located via 
the 'drop out' technique with recent UDF estimates from Q) and (0), (0), uncor- 
rected for extinction, as presented in (6). The integrated abundance to z ~10, while 
uncertain, appears to be insuffient for cosmic reionization. Possible explanations 
include an abundant population of lower luminosity sources or vagaries of cosmic 
variance in the small UDF field. Our surveys (§4 and §5) test both hypotheses. 



(SEDs) can be constructed and compared with population synthesis models 
to constrain the age, star formation history, and stellar mass of galaxies at 
z~6-7. In providing the opportunity to characterize the stellar mass density 
at z~6-7, Spitzer allows the first constraints to be placed on the star formation 
rate density integrated over the first 800 Myr of the Universe. 



The first demonstration of an established stellar population at high redshift 
was via the strongly-lensed z~6.8 galaxy in Abell 2218 (llh. The initial iden- 
tification of this triply-lensed source was based on detections in the F814W, 
F850LP, and F160W HST filters (0.801 /an, 0.905 /an, and 1.612 /an re- 
spectively). Subsequent detections were made at 3.6 /im and 4.5 /an with 
Spitzer for one of the components of the lensed system and at 1.1 /an with 
HST for two of the three lensed images. With the new photometric data, the 
galaxy's SED was constructed and compared to model SEDs generated with 
the GALAXEV population synthesis code (|14l). Best-fit SEDs of this source 
(Fig 2) are obtained by varying the galaxy age and redshift for a given star 
formation history and minimizing % 2 . The results suggest that the age of the 
galaxy is between 40-450 Myr, corresponding to a formation redshift of Zf ~ 
7 - 12. In each case the best-fit age of the galaxy is greater than the assumed 
e-folding star formation timescale, indicating that the galaxy would have been 
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Fig. 2. SED of the z ~6.8 lensed source in Figure 2 from (fllh . The combination 
of Spitzer and HST photometry spans the rest-frame UV to optical determining 
the luminosity-weighted stellar age (~100 Myr) suggesting a credible source of UV 
photons since z ~10. 

more luminous in the ultraviolet at earlier times during its active star-forming 
phase. Given the small survey area required for its discovery, such sources 
may have a reasonably high surface density suggesting that a population of 
intense star-forming galaxies (10-100 M yr _1 ) may be very common at z>7, 
contributing significantly to the UV photon budget required for reionization. 



The key feature which makes the lensed source in Abell 2218 of great diag- 
nostic value is the Balmer break observed between the NICMOS F160W and 
3.6 jttm bands; this indicates an age > 100 Myr and a declining star formation 
rate. Similar signatures have been subsequently seen in more massive systems 
at z~6-7 (Msteiiar ^ 10 10 M ). The SEDs of four i-drops in GOODS-South, 
spectroscopically confirmed by us at Keck, are constructed in (Jlif ). Two of 
the four galaxies have robust Spitzer detections at 3.6 fim and 4.5 /im. Com- 
parisons to population synthesis models (J14J) yield stellar mass estimates of 
2-4 x 10 10 M and ages between 250-650 Myr (z f ~ 7.5-13.5). Very recently, 
similarly massive systems have been presented for a larger sample in the UDF 
(J24l ). In order to build up such a stellar mass, the average SFR at earlier times 
must be greater than that observed at z ~6. 

Perhaps the most intriguing high-redshift object studied with Spitzer is a z~7 
gala xy in the Hubble UDF with a stellar mass claimed to exceed 10 11 M 
(|20j). The high-redshift interpretation is favored as a result of the strong break 
in flux between the K-band and 3.6 fim (interpreted as the rest-frame 3646 
A Balmer break) and the null detection in deepest optical broadband images 
ever obtained (3a detection limit of 31.5 in AB mags over a 1.5 arcsec aperture 
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in the combined BViz ACS images of the UDF). The best-fitting population 
synthesis models have very little dust extinction and ages of several hundred 
Myr, placing the epoch of formation of the bulk of this galaxy's stellar mass 
at Zf > 10. Even if star formation occured in several sub-galactic units which 
subsequently merged, remarkably high star formation rates (> 1000 M yr _1 ) 
may be required to assemble the stellar mass of this system by z~7. 



4 A Systematic Survey for Strongly-Lensed Lya emitters at z>8 

4-1 Survey Strategy 

The above discussion provides a clear motivation for searching for earlier ac- 
tivity and, prior to JWST and TMT, lensing offers a reasonable prospect of 
making some progress. We have undertaken a Keck spectroscopic survey for 
strongly-lensed star-forming sources at z>8. Using NIRSPEC in the J-band 
(sensitive to Lya emitters at 8.5<z<10.2), we scan the critical lines of ~10 
lensing clusters. By only observing along regions of very high magnification 
(x30 across most of the slit), we become sensitive to feeble sources produc- 
ing stars with rates as low as ~0.1 M yr _1 . All clusters used in our survey 
have spectroscopically-confirmed redshifts for a number of multiple images 
and arcs, allowing for the derivation of detailed magnification maps for high- 
redshift sources (Figure 3). 

Redshift identification of candidates is aided by several factors. The cluster 
lensing model constrains the redshift of any strongly-lensed source for which 
multiple images are detected via the observed image flux ratios and positions. 
In addition, by using NIRSPEC to search for Lya emitters in the same clusters 
studied with LRIS on Keck in the optical (0), we are able to correlate emission 
lines seen in the J band with existing optical spectroscopic data to clarify 
possible interlopers which are not multiply-imaged. Further, deep HST optical 
broadband imaging is available for all of our clusters. Lya emitters at z>8 
should be virtually invisible in the optical due to IGM absorption. Hence, the 
presence of optical flux spatially coincident with a candidate z>8 Lyo; emitter 
would indicate that the candidate is at low- redshift. 



4-2 Results and Discussion 

Twelve Keck nights have been allocated for this project. Nine nights have been 
completed thus far, of which 5.5 have been clear. During these 5.5 nights, 31 
slit positions were mapped toward seven clusters with an average integration 
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Fig. 3. The Keck NIRSPEC blind survey: slit positions observed toward Abell 68. 
Amplification contours (x5,30) are overlayed on the HST WFPC2 R-band image. 
The magnification appropriate for z>6 background sources is >30 for the entire 
extent of each of the three slit positions. 

time of 1.4 hours each. Each slit position consists of 6 to 9 ten minute exposures 
that were dithered parallel to the length of the slit. 

Spectra are flat-fielded and sky-subtracted with IDL routines following the 
optimal spectroscopic reduction techniques presented in (|2lh. We determine 
the pointing of each exposure on the sky by registering images from the NIR- 
SPEC slit-viewing camera to HST images of the same field. Offsets between 
exposures of a given slit position are calculated from the registered slit view- 
ing camera images, and the spectra are subsequently shifted and combined. 
In Figure 4, we present the reduced final two-dimensional spectra of three slit 
position toward Abell 68. 

The sky-subtracted two-dimensional spectra were independently inspected for 
candidates by the authors. A catalog of fifteen possible high-z Lya emitters 
was constructed and cross-checked with deep HST optical and ground-based 
near-IR images for corresponding broadband emission. Most candidates lack 
convincing emission in the ground-based near-IR broadband images. Final 
confirmation of these sources awaits deeper HST near-IR imaging (see § 5). 

Our limiting Lya line flux is defined to be the signal in an aperture 1.0 arcsec 
by 12 A that is five times the root-mean-square fluctuation in an aperture of 
that size. The spatial dimension was chosen to be roughly twice the typical 
seeing disk of the observations, and the spectral dimension was chosen to 
match the typical Lya line width of high-redshift galaxies, ~ 300 km s _1 
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Fig. 4. Left: Reduced two-dimensional spectra of three slit positions toward Abell 
68. A faint continuum trace from a foreground spiral galaxy is apparent in each of 
the three slit positions. Right: The minimum detectable flux (5a) for an observed 
Lya line as a function of wavelength (and redshift) from a 1.5 hour exposure with 
NIRSPEC. Depending on the local magnification, this observed limit of <10 -17 erg 
cm -2 s -1 corresponds to a star formation rate of ~ 0.03-0.1 Mq yr -1 . 

(0). Typical ha limiting line fluxes between OH lines range between / a =3 x 
10~ 18 erg cm" 2 s -1 and lx 10 -17 erg cm~ 2 s -1 (Figure 4). The corresponding 
characteristic limiting Lya line luminosity of our survey is ~10 41 erg s -1 after 
correction for the typical x30 magnification. Converting this line luminosity 
to a star formation rate, assuming 1 Mp yr _1 of star formation produces 10 42 
erg s _1 of luminosity in the Lya line ((22), we find the typical unlensed limiting 
star formation rate is ~0.1 M Q yr _1 . 

To obtain a rough estimate of the surface density of galaxies we would expect 
to find in our survey if low-luminosity galaxies dominate the reionization pro- 
cess, we consider two scenarios: 1) an extended reionization period (z=20 to 
z=7, At=575 Myr) with a moderate ionizing photon escape fraction (f c =0.5) 
and 2) a relatively quick reionization period (z=10 to z=7, At=300 Myr) with 
a lower ionizing photon escape fraction (f c =0.1). With these assumptions, we 
find that the surface density of sources in our survey must exceed 200 arcmin -2 
for the extended reionization case and 2000 arcmin -2 for the abrupt reioniza- 
tion scenario (Figure 5). Given our expected survey area of 0.015 arcmin 2 
(unlensed), upon completion of our survey we should detect several sources 
under a range of reionization scenarios. 
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We note that realistically, galaxies of all luminosities contribute to the reion- 
ization process. If the number density of galaxies continues to rise steeply 
toward star formation rates below the limiting rate of our survey, the num- 
ber density of star-forming sources required for reionization is lower than we 
estimate above. We consider this possiblity in more detail in a future paper 
(Stark et al. 2005, in prep). 

At the time of writing about half of our survey data has been analyzed sug- 
gesting that the abundance of galaxies at these early times is not significantly 
greater than that required for reionization under the assumptions mentioned 
above; however, it is clear from Figure 5 that our ability to detect reionization 
sources and provide meaningful constraints on the contribution of intrinsi- 
cally faint galaxies to reionization at z>8 will be significantly enhanced upon 
completion of our survey. 

Figure 5 also provides an opportunity for us to consider the challenge of ex- 
plaining the surface density of massive stellar systems being found by Spitzer 
at z ~6 0), (0), 0)- For each of the surveys (marked), we estimate the 



lower limit to the mass density from the sources so far seen and infer an equiv- 
alent lower limit to the (constant) star formation rate density required over 
the interval 6<z<20 to account for this mass. Not only are the estimates close 
to being adequate for reionization but, if the decline in Fig. 1 is correct, the 
presence of a few massive galaxies in the small fields so far surveyed already 
suggests the need for a very short intense period of activity at high redshift. 



5 An HST+Spitzer Search for Strongly-Lensed z=7-10 Galaxies 



To complement the Lya search, we have recently begun a HST large pro- 
gram (110 orbits) that will survey six cluster fields totaling 10 arcmin 2 in 
three wide-band filters (ACS-WFC /F850LP, NIC3/F110W, NIC3/F160W) to 
characterize the faint, star- forming galaxies responsible for the completion of 
cosmic reionization at z > 7. We achieve the necessary depth by boosting 
the sensitivity of HST with magnification from powerful lensing clusters at 
intermediate redshifts (Figure 6). 

We will select z > 7 sources using well-established color-selection techniques 
based on the neutral hydrogen break at rest-frame Lyman a (e.g., (Jig), (J5J); 

Incorporating the lensing magnification, more than half our survey area 
will reach unlensed star- formation rates of <1M yr _1 (for nominal source 
assumptions), i.e. better than the NICMOS UDF survey (J5j). Although lensing 
reduces the surveyed volume by the linear magnification, our total volume 
at UDF depths is still almost 25% of the NICMOS UDF field. Moreover, 
2.5 arcmin 2 is imaged more than one magnitude deeper in many sightlines and 
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Fig. 5. Estimates of the surface density of sources necessary for cosmic reionization 
over 7< z <10. The blue hatched area represents the likely range for various sce- 
narios; the lower bound represents an extended phase (7< z <20) and moderate 
ionizing photon escape fraction (f c =0.5) and the upper bound represents a rapid 
phase (7< z <10) and a low ionizing photon escape fraction (f c =0.1). Dashed curves 
illustrate our sensitivity limits (5a) to the surface density of star-forming sources 
using our present analyzed data (top curve) and that for the completed survey 
(bottom curve). A lower limit to the abundance of star- forming galaxies required to 
build the z ~6 massive galaxies described in (jlH ) .([lflh. |24h over the period 7<z<20 
is also displayed. Taken together with Fig. 1 these massive systems implies a very 
early period of intense activity sufficient for cosmic reionization. 

the fidelity of our sample will be assured by the well-understood geometrical 
effects of lensing. 

Both observational evidence (0) and theoretical conjecture (jjjj ) suggest the 
counts of faint z > 7 galaxies should be steeply rising at the UDF limit. 
Clearly at z ~ 7 only the most luminous galaxies are presently detectable (p). 
Extrapolating one magnitude deeper (Fig. 7) suggests we will reach a surface 
density of 11 arcmin 2 and find at least 5 sources fainter than the UDF limit, 
approaching the minimal luminosity density necessary for reionization. 
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Fig. 6. Deployment of HST NICMOS images (squares) for Abell 2219 - one of 
the target clusters. The outer curve defines the area within which, as a result of 
lensing magnification, the proposed survey will go deeper than the NICMOS UDF 
campaign. The regions inside the yellow curves represent areas of extreme (x 30) 
magnification. 

A major bonus of our program is that magnified sources are brought within 
reach of the Spitzer Space Telescope. All of our target clusters have been 
imaged by IRAC at 3.6 and 4.5 fj,m, to the same 5-cr point-source sensitiv- 
ity of 1 /iJj as Abell 2218. The extended wavelength baseline offered by the 
combination of HST and Spitzer can uniquely characterize physical properties 
(age, star formation history) for z~7 objects that would be too faint to oth- 
erwise be studied ffllh . By analogy with our analysis in Abell 2218, we can 
verify whether long-lived sources are the norm at z ~6-7. Tantalizingly, the 
combination of IRAC and NICMOS photometry also offers the prospect of 
searching for their more distant ancestors, z~8-10, via dropouts selected at 
the F110W/F160W boundary. 
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